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ABSTRACT 

Landing gear is one of the more prominent airframe noise sources. Techniques that diminish gear noise and 
suppress its radiation to the ground are highly desirable. Using a hybrid computational approach, this paper 
investigates the noise reduction potential of devices added to a simplified main landing gear model without small 
scale geometric details. The Ffowcs Williams and Hawkings equation is used to predict the noise at far-held 
observer locations from surface pressure data provided by unsteady CFD calculations. Because of the simplified 
nature of the model, most of the flow unsteadiness is restricted to low frequencies. The wheels, gear boxes, and 
oleo appear to be the primary sources of unsteadiness at these frequencies. The addition of fairings around the gear 
boxes and wheels, and the attachment of a splitter plate on the downstream side of the oleo significantly reduces 
the noise over a wide range of frequencies, but a dramatic increase in noise is observed at one frequency. The 
increased flow velocities, a consequence of the more streamlined bodies, appear to generate extra unsteadiness 
around other parts giving rise to the additional noise. Nonetheless, the calculations demonstrate the capability of 
the devices to improve overall landing gear noise. 

1 INTRODUCTION 

The past thirty years have seen significant reductions in jet noise through the adoption of high-bypass-ratio 
turbofan engines on civil aviation transports. Formerly unimportant noise sources such as the airframe have 
now become a major concern for noise certification and environmental considerations. Airframe noise is most 
important during aircraft approach and landing, when engines are operating at reduced thrust with high-lift 
devices and landing gear deployed. Wind tunnel tests 1-3 and fly-over measurements 4 have revealed the leading- 
edge slats, flap edges, and the landing gear to be the major contributors to airframe noise. Each of the three 
primary sources of airframe noise are important on different classes of airplanes, but the main landing gear 
is a dominant source on most modern wide -body transports. Although flow computations of high-lift devices 
such as flaps and slats received considerable attention in the last decade, the intricacies of a landing gear flow 
field and its associated sound sources have remained virtually unknown due to overwhelming geometrical 
complexities. Nonetheless, computational studies of landing gear are beginning to be performed . 5 Souliez et 
al . 6 used an unstructured grid technique to investigate the same baseline landing gear used in this work, but 
the calculations were restricted to laminar flow with the door removed. The authors previously investigated 7 
the baseline gear using a structured grid and the same numerical algorithm employed in the current work. 
In this paper, we present a computational analysis of noise reduction devices applied to the baseline gear. 
Our computational approach involves a hybrid strategy. In the first step we perform an unsteady Reynolds 
Averaged Navier-Stokes (URANS) simulation to provide a near-field solution. Despite continued advances 
in computational resources and numerical algorithms, it is still prohibitively expensive and often infeasible to 
attempt to resolve wave propagation from near-field sources to far-field observers. Integral techniques that can 
predict the far- field signal based solely on near-field input arc a means to overcome this difficulty. Hence, the 
Ffowcs Williams-Hawkings (FW-H) equation 8 solver described by Lockard 9 is used to predict the acoustic 
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signature at various observer locations using surface pressure data from the computational fluid dynamics 
(CFD) calculation. Additional details can be found in reference 10. 


2 ACOUSTIC EQUATIONS 


The FW-H equation can be written in differential form 1 1 as 


where 
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The dipole term F t involves an unsteady force, and Q gives rise to a monopole-type contribution that can 
be thought of as an unsteady mass addition. The function / = 0 defines the surface outside of which the 
solution is desired. The normalization |V/| = 1 is used for /. The total density and pressure arc given by 
p and p, respectively. The fluid velocities are u t , while the v r represent the velocities of the surface /. The 
Kronecker delta, dij, is unity for i = j and zero otherwise. The ambient speed of sound is denoted by c G . A 
prime is used to denote a perturbation quantity relative to the free-stream conditions denoted by the subscript 
o. The Cartesian coordinates and time arc x t and t, respectively. The usual convention, which is followed 
here, involves a quiescent ambient state with / prescribed as a function of time so that it always surrounds a 
moving source region of interest. H(f) is the Heaviside function which is unity for / > 0 and zero for / < 0. 
The derivative of the Heaviside function H'(f) = 6(f) is the Dirac delta function, which is zero for f =6 0, 
but yields a finite value when integrated over a region including / = 0. The inviscid paid, Pij = p6ij, of the 
compressive stress tensor P \j is used in this work. 

The FW-H equation is an exact rearrangement of the Navier-Stokes equations that allows one to determine 
the acoustic signal at distant observer locations if the details of the source region arc already known. Hence, 
the Navier-Stokes equations still need to be solved, but only where nonlinear and viscous effects arc impor- 
tant. All of the lineal - propagation can be determined by the FW-H equation. For three-dimensional flows, 
the time-domain FW-H formulations developed by Farassat 12 are efficient and amenable to numerical compu- 
tations. Some additional efficiency can be obtained by restricting the source to uniform, rectilineal' motion. 
Furthermore, the equation can be solved in the frequency domain which can be useful if one is only interested 
in analyzing certain frequencies. The frequency domain solution of FW-H equation can be written in the form 9 
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The quadrupole term is denoted by t Q and includes effects such as nonlinear propagation and refraction. 
In this work, the quadrupole contribution is expected to be small and is neglected. Souliez et alf performed 
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FW-H predictions of landing gear noise using solid and porous integration surfaces and found the solutions 
to be nearly identical in the far-held, although discrepancies were noted in the near held. The porous surface 
enclosed a signihcant region around the gear and should have included most of the quadrupole effects. In this 
work, our attention is restricted to the far-held where it should be valid to neglect the quadrupole. 

3 SIMPLIFIED LANDING GEAR MODEL 



(a) Baseline (b) Modihed 


Figure 1: Surface grids for a landing gear 



Figure 2: Details of landing-gear modifications 

The simulated geometry is a four-wheel landing gear model that approximately represents a ten percent 
scale Boeing 757 main landing gear. The model geometry is fairly complex and composed of four wheels, 
two side struts, an oleo, a side-door, yokes, a pin, and other structures that join the system together (hgure 1 
(a)). The gear assembly is mounted on a hat plate that represents the aircraft wing. The baseline structured 
grid consists of 155 blocks possessing a total of 1.8 million grid points. The modihed landing gear in hgure 
1(b) has 217 blocks and 2.1 million grid points. The modihed gear includes fairings around the gear boxes 
and between the wheels as well as a splitter plate behind the oleo. Figure 2 presents two views of the modi- 
hed landing gear with some of the components removed so that the added fairings are easily identihed. The 
purpose of the fairings is to smooth the how around the gear boxes and wheels in order to minimize bluff- 
body separation zones that are susceptible to large-scale how oscillations. The splitter plate 13 is supposed to 
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prevent vortex shedding from the oleo in the region where the door doesn’t already disrupt shedding. Figure 
1 also shows the grid distribution on the surface of the landing gears. The reference length scale is the gear 
wheel diameter (3.7in/0.09398 m), and the freestream Mach number is 0.2. The Reynolds number based on 
the wheel diameter is 1.23 x 10 6 . The CFD calculation employs the three-dimensional, time-dependent code 
CFL3D, 14,15 developed at NASA Langley Research Center to solve the three-dimensional, time-dependent, 
thin-layer Reynolds-Averaged Navier-Stokes (RANS) equations. A more detailed discussion of the CFD cal- 
culation can be found in the paper by Li et al . 16 The unsteady calculations were run extensively to eliminate 
all transients. The dimensionless time step based on the wheel diameter was c 0 At/l = 0.01. 




(a) Baseline 


(b) Modified 


Figure 3: Instantaneous perturbation pressure on landing gear. 

Contours of the instantaneous perturbation pressure fluctuations on the gear solid surfaces are displayed 
in figure 3. The surface pressure shows the footprint of the highly nonlinear and complicated interactive near- 
field flow dynamics. The pressure is nondimensionalized by p 0 Cg. Clearly, both the baseline and modified 
gears have significant pressure fluctuations, but the frequency and peak intensities cannot be inferred from 
instantaneous quantities. However, the high frequency oscillations on the door of the baseline gear are clearly 
absent for the modified gear. These waves emanate from small cavities between the door and the upper and 
lower yokes on the oleo. For the modified gear, these cavities were filled. 

4 NOISE CALCULATIONS 

The noise calculation for the baseline gear involves 181 total subsurfaces comprising the data surface. The 
modified gear has 224 subsurfaces. All of the subsurfaces are impenetrable, so only the pressure histories are 
needed. The subsurfaces are a natural consequence of the block structured grid used for the CFD calculation. 
Each subsurface is a boundary of one of the blocks comprising the grid. Because the problem is so large, the 
FW-H calculation would have to be performed on subsurfaces regardless of the grid topology. The complete 
time history record for all subsurfaces requires approximately 14 GB of disk space. Over 12,000 nondimen- 
sional time samples with c a At/l = 0.02 have been collected. Five segments of 4096 samples are used to 
perform the noise calculations. The results presented in the paper represent an average over the five segments. 
The segment length corresponds to the time for the flow to pass a wheel sixteen times. 

A Ffowcs Williams-Hawkings solver 9 written specifically for airframe noise applications is being used to 
perform the noise calculations. The calculated sound pressure level spectrum for an observer located directly 
underneath the gear is presented in figure 4. The observer is located 100 wheel diameters directly beneath 
the gear. Most of the noise is concentrated at low frequencies because of the absence of the smaller scale 
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Figure 4: Spectra of the pressure for an observer 100 wheel diameters directly below the gear. 


subcomponents in the model and possible excessive diffusion caused by the turbulence model or insufficient 
grid resolution. However, two high-frequency tones are evident in the signal as shown in figure 4. These tones 
are those caused by the resonances in the cavities between the yokes and the door. The figure shows that the 
tone from the upstream cavity is around 20 kHz, and the downstream cavity generates a tone at 26 kHz. There 
is no evidence of the tones in the modified gear spectra. 


4.1 SUBSURFACE NOISE PREDICTIONS 


Although much of the smaller scale detail is missing from this gear model, the frequency content below 6 
kHz is relatively rich and it is not apparent which components arc contributing to the different portions of the 
spectrum. There are several distinct tones that are probably the result of shedding from different components. 
To investigate the source of the different peaks, the landing gear was divided into 1 1 regions as show in figure 
5. Each region is colored differently to identify each of the subdomains. 
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(a) Baseline 


(b) Modified 


Figure 5 : Landing gear colored by component. 


The prediction using all of the surfaces (Gear + Ceiling) is compared with the results when using only the 
ceiling in figure 6. The results are for an observer located 100 wheel diameters directly below the gear. The 
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results do vary with observer position, but the general trends arc similar at most observer locations of interest. 
The figure shows that for both gears the ceiling accounts for almost all of the noise. The signal from the 
components is 5-10 dB lower than that of the ceiling. Although one would expect the ceiling to be important 
because of reflections, it should not be dominant. A calculation without the ceiling where all of the gear 
components are mirrored still does not produce levels anywhere near those observed with the ceiling alone. 
Clearly, the ceiling is acting as more than just a reflector as is borne out by figure 7 that shows the instantaneous 
perturbation pressure contours on the gear and ceiling looking from below. The most intense fluctuations occur 
on the ceiling in the wake of the sidebars. Although it is likely that the wake from the sidebars would interact 
with the ceiling, the fluctuations are actually being amplified across and along block boundaries because of 
insufficient grid resolution. Each of the discontinuities in the contours in figure 7 represents a block interface 
where extensive patching is used to prevent the propagation of fine resolution to unimportant regions and 
thus reduce the number of grid points. Because it cannot be determined how much unrealistic amplification 
is occurring, the prediction when the ceiling is excluded will be used as a reference in the remainder of the 
paper. The spectra comparison using just the gear components is shown in figure 8. Clearly, the landing gear 
modifications have produced significant broadband noise reduction, but two additional tones arc also apparent. 
Furthermore, the impact of the modifications on the peak of the spectra has been minimal. A subsurface 
analysis and source identification will be used to help identify the causes of the tones and the noise reduction. 
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(a) Gear + Ceiling 


(b) Ceiling 


Figure 6: FW-H results for an observer 100 wheel diameters directly below the landing gear in a M = 0.2 
flow. 


The predicted noise for each of the subsurfaces in figure 5 is shown in figures 9 and 10. The peak levels 
for the 4.4 kHz tone and its harmonics appear to be associated with the oleo (figure 9(f)). The tones appear in 
the spectra for the other components because the signal is reflected off the other parts. The tone is associated 
with unsteadiness in the region around the yoke on the oleo. The tone is also apparent in the baseline spectra 
at 3.7 kHz. Apparently, streamlining the model accelerated the flow velocity around the oleo increasing the 
intensity of the unsteadiness and slightly changing the frequency. Although one can detect vortex shedding 
in this region from the CFD solution, it is asymmetric and associated with a complex interaction of the flow 
around the oleo, yoke, and door. The splitter plate appears effective at eliminating normal Strouhal shedding 
which is expected around 1 kHz for the contracted region of the oleo where the flow is not disrupted by the 
door. The 1 kHz tone observed for the baseline model is virtually eliminated by the splitter plate. 

The broadband shape of the gear spectra is generated by the gear boxes (9(c)), the connectors between the 
gear boxes (9(d)), and the wheels (10(c)). Significant broadband reduction is apparent for all three components 
on the modified gear. However, the reduction below 500 Hz is negligible, so the overall peak of the noise is 
unchanged. 
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M=0.2 


(a) Baseline 



M=0 . 2 


(b) Modified 


Figure 7: Instantaneous perturbation pressure on the landing gear. 



Figure 8: FW-H results for an observer 100 wheel diameters directly below the landing gear in a M = 0.2 
flow. Only the gear components are used to make the prediction. 
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Noise reduction is also observable in the spectra for the axles (9(a)), door (9(b)), pin (10(a)), and yokes 
(10(d)). In some cases, the observed reduction is primarily a result of reducing sources on other components 
whose footprint extends to other areas. Nonetheless, the general trend of reduced noise is encouraging. The 
spectra for the sidebars (10(b)) is relatively unchanged by any of the modifications which is expected because 
they arc unchanged from the baseline configuration. However, the spectra for the hub caps (9(e)) actually 
increased, presumably because of accelerated flow around the wheels caused by the additional streamlined 
parts. 

5 RESOLUTION AND TURBULENCE MODEL 

The spectra obtained from the simulations is restricted to low frequencies which makes it difficult to assess 
the total effect of the noise reduction devices. The model is approximately 10% of full scale, so full scale 
frequencies up to 600 Hz have been obtained from the simulations. Measurements of landing gear noise 
typically extend to several kHz. One possible cause of the loss of the high frequency content is the mesh. It 
is difficult to maintain adequate resolution throughout the flowfield while still making the problem tractable. 
Furthermore, excessive patching is required for a complex geometry such as a landing gear, and errors at block 
interfaces arc apparent in the solution. However, a calculation on a finer grid has failed to completely fill 
out the spectra. Another possible source of error is the turbulence model. We have used an unsteady RANS 
calculation with a k — u> turbulence model tuned for steady flow. Although qualitatively good results have 
been obtained with this approach for unsteady problems, it is well known that the turbulence model is overly 
dissipative. We have also run the landing gear problem as a detached-eddy simulation 17 (DES) as proposed 
by Spalart. The DES model essentially reduces the level of eddy viscosity in regions away from solid surfaces 
when the grid is sufficiently fine. The comparison of the spectra in figure 1 1 shows that the DES solution 
does have frequency content out to 9 kHz, and the levels are significantly higher above 2 kHz. Although a 
significant improvement, the levels arc still lower than expected based on experiment. 

6 CONCLUSIONS 

The current investigation has validated the noise reduction potential of fairings around the wheels and gear 
boxes of a landing gear, as well as the potential elimination of vortex shedding behind cylindrical bodies by the 
addition of a splitter plate in the wake. Although significant broadband noise reduction was obtained, the peak 
of the spectrum was unchanged. Furthermore, the level of a tone was significantly increased. The tone was 
affected because the flow around certain components was accelerated by the streamlining of some of the bluff 
bodies comprising the gear. Hence, one should be aware that there can be unexpected consequences of the 
addition of noise reduction devices. Two high frequency tones associated with unsteadiness in small cavities 
between the yokes and the door were eliminated by filling the cavities. Small cavities should be avoided when 
possible because of the potential for resonances in these regions. 

The accuracy of the simulations was examined through grid refinement and changing the turbulence model. 
A finer mesh only resulted in minor changes to the calculation, but a DES simulation did improve the range 
of resolved frequencies. However, there is still a large discrepancy with what is observed in experiments. 
Nonetheless, it is expected that the effect of noise reduction devices investigated in this work arc most influ- 
ential in the range of frequencies resolved by the current calculations, so the most important aspects of the 
physics have been captured. 
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Figure 9: Noise contributions for an observer 100 wheel diameters directly below the landing gear in a M 
0.2 flow. 
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(a) Pin (b) Sidebars 




(c) Wheels (d) Yokes 

Figure 10: Noise contributions for an observer 100 wheel diameters directly below the landing gear in a 
M = 0.2 flow. 



Figure 11: Spectra comparison for an observer located 100 wheel diameters directly below the gear in a 
M = 0.2 flow. 
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